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Abstract 

The one-loop Wilsonian renormalization group flows of the two dimensional two-flavors Gross- 
Neveu model with minimal doubling fermion are calculated numerically. The off-diagonal mass com- 
^ . ponents which are non-Hermiticity are generated by the renormalization group flows. We consider 
a relation among 75-Hermiticity, R-Hermiticity and PT symmetry. We point out that satisfied two 
of the three conditions is a sufficient condition for satisfied another condition but not a necessary 
condition, this relation has a connection with Hermiticity. In the case of periodic and continuum 
kinetic term, both 75-Hermiticity and R-Hermiticity should be satisfied for consistency. 

Because kinetic terms for the minimal doubling fermion do not have R-Hermiticity, non-Hermiticity 
effective masses appear. 
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1 Introduction 



Lattice gauge theory is one of powerful tools to reveal non-perturbative quark dynamics [lj. In 
this formulation, space-time coordinate is discretized and physical variables are defined at sites and 
links. We can calculate observables in strong coupling region using some techniques, for example 
high temperature expansion. Especially Monte-Carlo simulations are effective methods to investigate 
non-perturbative physics and carried out actively. 

As well known, the doubling problem that fermion has redundant degrees of freedom called "dou- 
blers" is a main problem to apply lattice formulation to theory involved matter fields, for example 
quantum color dynamics (chromodynamics) (QCD). It is impossible to remove doublers unless we give 
up some symmetries or properties because of the no-go theorem of Nilsen and Ninomiya [3J . To over- 
come this problem, many fermion formulations are constructed, for example Wilson fermion [1] and 
KS fermion [2j. The exact chiral symmetric action is important symmetry to analyze non-perturbative 
QCD, however this symmetry is incompatible with removing doublers. 

Recent year, Creutz constructed exact chiral symmetric lattice fermion action [4], and Borici fitted 
it to orthogonal lattice [5]. Few decades ago Karsten constructed a fermion formulation which has 
same structure and is different form from Creutz action [6j. These fermions are called "minimal 
doubling fermion" . The minimal doubling fermion actions break (hyper-)cubic symmetry and some 
discrete symmetries, such as charge conjugation(C), parity transformation (P), time reflection(T) and 
so on. Many properties of the fermions are analyzed [Bj-PS] 0. 

In this paper, we apply minimal doubling fermions to the two dimensional N-flavors Gross-Neveu 
model which is a toy model of QCD. We calculate the one-loop Wilsonian renormalization group flows 
(RGFs) in this model. 

As a result, the off-diagonal mass components which are non-Hermiticity are generated in this 
flows, even then starting from massless initial conditions. From this simple analysis, we point out 
that minimal doubling fermions have a serious defect, non-Hermiticity. Then we discuss the relation 
among R-Hermiticity, which is defined in Section 4, PT (or C) symmetry and 75Hermiticity. 

And we also suggest a criterion to forbid this phenomenon. 

This paper is constructed as follows. In Section [2] we will describe for minimal doubling fermion 
action in two dimension. In Section [3] we will briefly review N-flavors Gross-Neveu model, and then 
numerically analyze one-loop Wilsonian renormalization group flows around trivial fixed point using 
naive action and minimal doubling actions. In Section [4] we will discuss Hermiticiy, PT symmetry 
and 75Hermiticity conditions for a general lattice fermion kinetic term. In Section [5] we will conclude 
and discuss about this relation and minimal doubling fermion. 

2 Minimal Doubling Fermion Action 

In this section, we will review minimal doubling fermion actions briefly [4]-|17|. 
We define kinetic terms of naive action(NA) and two minimal doubling actions(MDAs) in two 
dimensional momentum space as follows: 





1 See ref. |16lfT7] about properties and symmetries of hyperdiamond lattice minimal doubling fermion. 
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where the subscript "kin" means a kinetic term, and 
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(2.2) 



We fix a value of lattice space, a = 1, from now on. The subscripts 1,4 means space and time 
components respectively. The -D m di(j?) and -D m d2 are called "twisted ordering action" and "dropped 
twisted ordering action" respectively [13J. In two dimension, the NA has four zero-modes and the 
MDAs have two ones which appear in the following momenta: 



Andl 
D m d2 



p= (0,0), (0,vr), (vr,0)and(^,vr), 

p= (0,0)and(vr/2,7r/2), (2.3) 
p= (0,0) and (0,vr). 

D(p + q) with D{p) = 0, and they 



The "doublers" appear around each zero- modes as Dip) 
contribute observables. 

In the case of NA, the half doublers that half doublers have same chirality and the others have 
opposite. 

In cases of MDAs, they have opposite chirality each other. The NA and MDAs have 75Hermiticity: 



75^>(p)75 = D\p). 
For massless, they also have chiral symmetry: 

j 5 Dip) + Dip) l5 = 0. 



(2.4) 
(2.5) 



The MDAs violate (hyper-)cubic symmetry and some discrete symmetries. We define charge conju- 
gation^), parity transformationfP), time reflaction(T) and their combinational transformation laws 
acting on a fermion kinetic terrro 
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We present these symmetric properties of the NA and MDAs in Tab ]Q 



3 N-flavors Gross-Neveu model and Renormalization Group Flow 
in Two Dimension 

In this section, we describe the N-flavors Gross-Neveu(GN) model |20] and calculate Wilsionian 
renormalization group flows(RGFs) using the NA and MDAs numerically. Firstly we will review the 
N-flavors GN model and then we will calculate the RGFs. 

2 We can apply same laws to four dimensional theory, replacing with pi —5- p. 

3 "T" represents site and link reflection in lattice space. In the case of the -D m d2, it has link reflection positive but 
not site reflection positive [8] . 
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e 1: Discrete symmetry for the NA and the MDAs 
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3.1 Action and Symmetry 

We define the continuum Euclidian Lagrangian of N-flavors GN model as follows: 

2 

£ G n = $(d-7 + m)^- — (4>4>) 2 , (3.1) 

where m is a fermion mass and g is a coupling constant of four fermi interaction. We omit flavor 
indices if we do not have to write explicitly, ipip = ^iLi'ipiipi, where "i" means flavor degrees of 
freedom. 

This Lagrangian has U(l) symmetry: 

V> ->• e i( V, 

i> ^e~ ie . (3.2) 
In the case of massless fermions, this Lagrangian has chiral Z4 symmetry: 

i> -)• VKnsf. (n = 0,l,2,3) (3.3) 

In the case of massive fermions, chiral Z4 symmetry reduces to chiral Z2 symmetry (n = 0,2). In 
addition, if all flavors have same masses it has the SU(N)p symmetry: 

tpi ->■ Uijipj, 

A -> jjUb (UeSU(N)) (3.4) 
It is convenient to redefine the GN action using an auxiliary scalar field a instead of (ipip): 

>Cqn = i> (d ■ 7 + m) ip + ycr 2 + gaipi/j. (3.5) 

According to this manipulation, we can obtain the action which involves Yukawa interaction instead of 
four fermi interaction. According to perterbative calculation, the NG model has asymptotic freedom 
property [2T1I22]. 

3.2 Wilsonian Renormalization Group Flow 

Using the Wilsonian method |23j . we calculate numerically the RGFs for the mass and coupling 
constant starting from the trivial fixed point (m = g 2 = 0). 

We explain how to calculate Wilsonian RGFs in Appendix lAl 

In the case of MDAs, we use each zero-modes as the different flavor fermions, and in the case of 
NA we use only two zero-modes, p = (0, 0) and (it, it). We now represent the spinor indices explicitly, 
and we distinguish 0, 1 from 2, 3 as different flavors. We assume that high frequency modes of fields 
are not effective, therefore we neglect their contributions, < \k\), -0(1 < |fe|) and a(l < \k\). 



We choose initial conditions for the mass as moo = m n = m 22 = m 33 = 0, ±0.25, ±0.5, and for 
the coupling constant as g 2 = 0, 0.2, 0.4. 

The off-diagonal mass components equal to zero in all cases. We will calculate numerically up to 
one-loop quantum effects and RGFs which toward IR from the initial conditions @. In our calculation 
we define 7-matrices as follows: 

i)>-»- (?;)• w> 

The results are referred in figHJ The RGFs of the NA and MDAs are similar forms and differences 
among each values are O(10 -3 ). 

In cases of MDAs, however, off-diagonal mass components are generated by the RGFs, except the 
initial value which is trivial fixed point(m = 0, g 2 = 0). The obtained results do not have Hermiticity. 

We show this fact in figj2]and fig{3] whose initial conditions are m = 0, g 2 = 0.2. The figEJa) and 
(b) are the RGFs of -D m di(p) and D m d2{p) respectively. The figGta) and (b) are relations between 
the off-diagonal mass components and iterations using -D m ai and -D m d2 respectively. Though the 
off-diagonal mass components amplify as the flows go to IR because of scaling effect, they do not 
break chiral Z4 symmetry. This fact means that the MDAs do not have Hermiticity and it seems to 
have any problems. In the following section, we will discuss this phenomenon and Hermiticity. 



4 75-Hermiticity, R-Hermiticity and PT symmetry 

In the previous section, we observed non-Hermiticity contributions in fermion effective mass com- 
ponents. In Euclidian formulation, a fermion kinetic term in continuum limit does not have Her- 
miticity but anti-Hermiticity. The NA also has an anti-Hermiticity kinetic term. In the case of NA, 
however, an effective mass does not have anti-Hermite components explicitly, therefore the cause of 
them is clearly added even function terms, (cosp — 1)7. We will consider Hermiticity and PT sym- 
metry, and what added terms generate explicit non-Hermite effective coupling constants. Hermicitity 
is discussed in ref.[18] at the case of the minimal doubling fermion on the Hyperdiamond lattice. We 
will focus on only kinetic terms from now on and discuss In the case of two dimension □. 

We will focus on only kinetic terms from now on. In the case of two dimension, we define a 
translation invariant kinetic term in momentum space as follows: 
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D(k) 



(4.1) 



(4.2) 



M =l,4 



where fu(k) are com 
and PT symmetry^ 



jlex numbers in general. We define three conditions, 75-Hermiticity, R-Hermiticity 



75 Hermiticity 
R-Hermiticiy 
PT symmetry 



D(k) = l5 D\k) 75 , 
D(k) = D\-k), 
D{k) = 75 .D(-%5- 



(4.3) 
(4.4) 
(4.5) 



4 According to estimate integrating part of one-loop calculation we used the sectional measurement method. We took 
the length of division Ap M = 0.01 (O(0. 01 2 ) error in the integrating part and it seems that accuracy is not very well. 
5 We can apply this discussion to the case of any even dimension. 

6 In following discussion, we can use C symmetry instead^ of PT symmetry because of CPT theorem. 
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Figure 1: The RGFs of MA and MDAs. Initial parameters are mu = 77122 = "^33 = ^44 = 
0, ±0.25, ±0.5, g 2 = (a)0, (6)0.2, (c)0.4. The RGFs toward IR which have large mass. 

These conditions are not independent of one another. We can easily lead to the fact that another con- 
dition is automatically satisfied if a kinetic term satisfy two of the three conditions. This is a sufficient 
condition but not a necessary condition. If f^{k) is pure imaginary, 75-Hermiticity guarantees anti- 
Hermite condition for a kinetic term. R-Hermiticity is also used as (anti-)Hermite condition, for ex- 
ample in ref.[19], however it is not well-defined because forward-derivative, -Dfd(fc) = (e* fc '' — l) 7^, 
satisfy this condition. 

We will show that R-Hermiticity is a condition for real effective coupling constants. We assume 
that a fermion kinetic term has R-Hermiticity and effective coupling constants is the following form: 

00 

Seff = go + J2 l(n) > ( 4 - 6 ) 

n=l 

with 

iin) = /n^-^..^(-^---.-^-n^(*i)' w 

where go is a real bare coupling constant, g e ff is an effective coupling constant, S a p{k) is a fermion 
propagator and 1^ is n-loop quantum effect which is constructed from r-fermion propagertors. We 
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Figure 2: The off-diagonal mass and coupling constant of the MDAs, (a)D m( ji, (b)-D m d2- Initial 
parameters are mn = 77122 = ^33 = "144 = 0, g 2 = 0.2. The RGFs toward IR which has large mass. 
The RGFs have irregular forms and the off-diagonal mass components are non-Hermiticity. 

act Hermite conjugate to the second term of r.h.s. The effective parameter is real if the following 
condition is satisfied, 

^k-arA.C-*!'-" >-M=j£i.^(*i," - ,M- (4-8) 

Ordinary Euclidian actions are constructed from Hermite terms except a fermion kinetic term, so 
that /( n ) is also constructed from only Hermite ingredients. Therefore we can use R-Hermiticiy as 
a reality condition for parameters, as long as (j4.8[) is satisfied. If 75-Hermiticity is not satisfied and 
R-Hermiticity is satisfied, one-loop fermion propagator produces diagonal mass components. In the 
continuum Euclidian action case, these components are not produced. Therefore both 75-Hermiticity 
and R-Hermiticity should be satisfied for consistency. 

Next, we will show that PT symmetry is always broken if we add extra kinetic terms to a NA to 
reduce to doublers 0. 

Statement 

In even dimension, a PT symmetric kinetic term which is assumed periodicity and continuity always 
has more than or equals to 2 d poles. 

Proof 

For simplicity, we also assume translation invariant H. A general 2ir periodic and continuum D{k) is 
the following form: 

d 00 

cos(n u k u ) + (C M (n) + iD^n)) sm(n u k v ) + E^] 7^, (4.9) 

where Au(n), B^(n), C^(n), D^(n) are real numbers and are complex constants. From PT sym- 
metry, 

A tl (n) = B ll (n) = E ll = 0, for all fi,n. (4.10) 



7 This argument has discussed in ref.[9] but not mathematically. 

8 We can similarly lead to same statement without translation invariance. In the case of non-translation invariant, a 
kinetic term is D(k,p) form and at least 4 d doublers appeaj;. 
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Figure 3: The off-diagonal mass and iteration of the MDAs, (a)-D m( n, (b)D mc i2- Initial parameters 
are m\x = m22 = "mm = 77144 = 0,g 2 = 0.2. The off-diagonal mass components are generated and 
they have non-Hermiticity. 

The D(k) always has two poles at k = and n for each dimension. 
Therefore D{k) has equal or more than 2 d poles. 

This statement equals to that we can not reduce to the number of doublers using PT symmet- 
ric kinetic terms [^|. In numerical simulation context, 75-Hermiticity is a very important condition 
to avoid sign problem. Assumed translation invariance, R-Hermiticity is not satisfied if D(k) sat- 
isfy 75Hermiticity and does not PT symmetry. Therefore effective parameters have explicit non- 
Hermiticity. In a process of rewriting from Minkovskian to Euclidian, Hermite fermion kinetic terms 
transmute anti-Hermite ones. A general lattice fermions which is reduced the number of doublers 
have possibilities of generating this anti-Hermite effective coupling constants. The R-Hermiticity is a 
criterion to remove non-Hermiticity. Because MDAs have only 75-Hermiticity, the fermion effective 
mass which is constructed from odd fermion propagators has non-Hermiticity components explicitly 
in RGFs. 



5 Conclusion and Discussion 

We have investigated the one- loop Wilsonian renormalization group flows (RGFs) of the two- 
flavors Gross-Neveu(GN) model in two dimension with naive action(NA) and two minimal doubling 
actions (MDAs). We observed that in cases of MDAs off-diagonal mass components which are non- 
Hermiticity are generated even massless initial conditions. In order to understand the reason for non- 
Hermiticity mass components, we considered 7sHermiticity, R-Hermiticity and PT (or C) symmetry 
conditions. These conditions are not independent, satisfied two of the three conditions is a sufficient 
condition for that another condition is satisfied automatically. However it is not a necessary condition. 
We suggested that R-Hermiticity is a condition for removing non-Hermicity. However, both 75- 
Hermiticity and R-Hermiticity should be satisfied. Becausethe diagonal mass components which is 
not consistent are generated if only R-Hermiticy is satisfied. Therefore we can not reduce the numbers 
of doubler with modified kinetic terms which have periodicity and continuity because of a relation 
between doublers and PT symmetry. 



9 We can not apply this statement to non-7 M linear case, for example Wilson fermion. 

8 



We applied this relation to the MDAs which have only 75-Hermiciticy. Because of the non-R- 
Hermite kinetic terms, explicit non-Hermite effective mass components appear. 

In lattice perturbative, theory we have to care whether we can remove relevant and marginal 
symmetry breaking terms with counter-terms 0. In massive two-flavors GN model, non-Hermite 
masses are small in comparision with diagonal masses and we can remove them with counter-terms. In 
massless case, non-Hermite masses appear however they do not break chiral Z4 symmetry. Therefore 
it is not so trouble at one-loop level. In more higher loop or non-perturbative case, we do not know 
how far parameters are involved the non-Hermite effects and should be fine-tuned. 
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This issues are discussed in the case of four dimensional gauge theory in ref.|ll] for example. 



Appendix 

A Wilsonian Renormalization Group 

In this appendix, we review a method to calculate the Wilsonian renormalization group flow in 
the case of GN model in two dimension [23] . 

We define the partition function of the GN model in momentum space q 



with 



where 
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and £gn is given in (|3,5p . We can treat of AT as a mass parameter of auxiliary field a. Here we assumed 
that the high frequency modes had already integrated and they do not contribute effectivelly. Then 
we split the field configurations as follows: 



where 



5 

4 

(?h(p) = cr(p) if 7 < bl < f' 
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(A.6) 



and the other fields are also split similarly [3 We choice renormalization conditions as follows: 
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where are renormalized i-point functions, itir, Nr, gR are renormalized parameters and arguments 
of r« are external momenta. In order to obtain effective parameters, we calculate one-loop effect 
and integrate out only high frequency modes: 
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1 We omit the subscript which means flavors. 



2 On account of numerical efficiency, we choose a division^ which split between ai and an as p = | 



where S(k) and D{k) are propagators of each fields presented below, "tr" is a trace operation of the 
fermionic indices and rj^p and rja- are rescaling parameters of fermion and auxiliary field respectively. 
We can define these parameters with dimensional analysis in the following values: 

/5\ 3 / 2 
5 

* = ? 

We can obtain propagators from the GN action: 



(A.15) 
(A.16) 

where Df{k) is one of the lattice fermion kinetic terms in (|2,2p and A; is a zero- mode momentum in 
(|2.4p . Substituting (|A.15|) and (|A.16|) to (|A.11[) - (|A.12|) . we can obtain effective mass and coupling 
constant after integrating out. 
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